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Hart et al. [9] 


rif — 2, partially quenched, 


1.0(2) 


Prelovsek et al. [12] 


n/ = 2, unquenched, 


1.58(34) 


Prelovsek et al. fl2] 


partially quenched 


1.51(19) 


Burch et al. [5] 
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Lattice QCD with Nf — 2 flavours of sea quark is used to explore the spectrum and decay of 
scalar mesons. We are able to determine the fei - ao mass difference and this leads to the conclusion 
that the lightest non-singlet scalar meson (ao) has a mass of 1.01(4) GeV. We determine from the 
lattice the coupling strength to KK and nrj. We compute the leptonic decay constant of the lightest 
non-singlet scalar meson. We discuss the impact of these lattice results on the interpretation of the 
Ho (980) state. We also discuss Kq states. 

PACS numbers: 12.38.Gc, 14.40.Cs, 13.25.-k 



I. INTRODUCTION 

The scalar mesons known experimentally do not fit into 
a tidy pattern, as found for vector or axial mesons, for 
example. Because the scalar mesons have S-wave decays 
to light two-body states (two pseudoscalar mesons), then 
the impact of these two-body channels on the scalar me- 
son can be sizeable. Thus the scalar mesons may have 
qqqq as well as qq components. For example, there are 
two ao mesons, ao(980) and ao(1450), known [j|. The 
ao(980) meson is closely associated with the KK thresh- 
old and it has been suggested that this is a molecular 
state. This is can be explored using lattice techniques. 
A further complication is that the flavour singlet scalar 
mesons can mix with scalar glueballs, although here we 
restrict our investigation to the flavour non-singlet scalar 
mesons from lattice QCD. 

There has been a long history of studying the scalar 
non-singlet mesons on the lattice. These states tend 
to have a poorer signal to noise ratio than the S-wave 
mesons 0] such as the p and tt, hence are less commonly 
studied. Much of the early literature on light P-wave 
mesons focussed on designing good interpolating opera- 
tors to create the mesons 0, Q ■ 

The quenched studies of the ap were complicated by 
the discovery of a ghost state that made the correlator 
for the ao particle, which should be positive definite in 
a unitary quantum field theory, go negative |j|. If this 
effect was not taken into account then the chiral extrap- 
olation of correlators was unreliable. Modern studies of 
this state such as those by Burch et al. [^ correct for 
the effect of the missing contribution to the 0++ corre- 
lator from the r]' meson. Prelovsek p\ has also studied 
the ghost state in the ao correlator using 2H-1 dynamical 
staggered fermions and mixed (chiral valence and stag- 
gered sea) fermions. In both cases, which are essentially 
partially quenched, deviant features are discovered. 

The non-singlet scalar mass is an input into the study 
of mixing with rfueballs in the singlet sector by Wein- 
garten and Lee Q. The ghost state was not taken into 
account and this led to problems with the chiral extrapo- 
lation of the non-singlet 0"'""'' meson masses. This mixing 



TABLE L Some results for the mass of the ao meson from 
quenched and partially quenched QCD that include the effect 
of the ghost state |J|. 



has also been discussed [a, |y| using unquenched lattices 
which avoids this problem. 

Alford and Jaffe |l3| used quenched QCD with ifq^ 
operators relevant to 0++ mesons. Their study claimed 
to see evidence for bound states in the q^q^ channel rel- 
evant to 0++ states. The work of Alford and Jaffe [13 
can be criticised for not taking into account the quenched 
ghost in the oq correlator. Only a subset of the correla- 
tors required for the singlet channel were computed. This 
is, perhaps, consistent in quenched QCD but clearly im- 
portant physics is omitted. 

The scalar collaboration are startin g to use lattice 
QCD techniques to study the k particle \U\ . 

Prelovsek et al. [ij extended the work of Bardeen et 
al. 4] on the effect of the ghost state in the ao channel 
to the partially quenched theory. By restricting lattice 
study to valence quarks heavier than the sea-quarks. Hart 
et al. Q were able to extrapolate to light quarks with 
no ghost contributions, obtaining an estimate for the ao 
mass of 1.0(2) GeV. 

In table m we collect together some recent numbers for 
the mass of the ao mass from some modern lattice cal- 
culations that take into account the ghost term. None of 
the calculations in table |l] had complete control over all 
systematic errors, such as finite size effects or the contin- 
uum limit, even within quenched QCD. The results for 
the lightest 0++ meson are mostly around 1.5 GeV. As 
we note above, the ao decays via the strong interaction, 
so a quenched QCD calculation may give a poor estimate 
of the particle mass. 



The MILC collaboration reported evidence for oq de- 
cay on the lattice in an unquenched lattice QCD calcu- 
lation with 2 + 1 flavours of improved staggered fermions 
with a lattice spacing of 0.12 fm 13]. In MILC's first pa- 
per they found the ag mass to be significantly lower than 
the mass of the bi and ai mesons. This was different 
behaviour from the quenched study with the same pa- 
rameters. MILC ;13] found that, for lighter quarks, the 
mass of the Oq meson was close to the sum of the tt and 
r] masses, where the mass of the rj was estimated using 
the Gell-Mann-Okubo formula. As the MILC collabora- 
tion |l4| ran unquenched calculations with even lighter 
sea quarks they confirmed that the lightest state in the ao 
channel lay below the irrj threshold. Using independent 
techniques on a subset of the configurations from MILC, 
Gregory et al. jl^j also found that the lightest state in 
the ao channel was below the ttt] threshold. Prelovsek j^ 
has studied oq decay using staggered chiral perturbation 
theory, concluding that taste violations in the staggered 
fermion formalism allow a small amplitude for the de- 
cay of the Oq state to two pions. The decay Oq —^ tttt is 
forbidden in the real world because of G parity. 

Since the current state of lattice investigation of scalar 
mesons is incomplete, more work is needed. In order to 
make a start in establishing the nature of scalar mesons 
from first principle in QCD, we address here the flavour 
non-singlet scalar mesons. As is well known, lattice QCD 
in the quenched approximation is not a consistent theory 
and this manifests itself as ghost contributions to the 
scalar meson propagation - arising from the spurious low- 
lying threshold in the ttt] two-body channel. We use here 
Nf — 2 dynamical gauge configurations so that we have a 
consistent field theory. The physical case, however, also 
has another light quark (the s quark) and has lighter u 
and d quarks than we are able to use on a lattice. Thus 
some extrapolation will be needed to obtain consequences 
for the physical spectrum. 

As we approach the limit of physical light quark 
masses, the scalar mesons become unstable: they are res- 
onances. On dynamical lattices these decay channels are 
open. Thus we need to have methods to cope with unsta- 
ble particles on a lattice. The study of hadronic decays 
from the lattice is not straightforward - see ref. [Tal . It is 
possible, however, to evaluate the appropriate hadronic 
matrix element from a lattice if the transition is approxi- 
mately on-shell. This allows us to estimate decay widths, 
provided that the underlying coupling is relatively insen- 
sitive to the quark masses. We follow methods generi- 
cally similar to those used by us to study p decay [l3| 
and hybrid meson decay [Ig. 

As well as the hadronic decay, one can also define a de- 
cay constant analogously to that defined for the weak de- 
cay of pseudoscalar mesons. We discuss the relevance of 
this and the determination from the lattice of the scalar 
decay constant. 



II. SPECTRUM 

As a by-product of our study of hybrid mesons, we 
have accurate lattice measurements of the ao, 6i and ai 
mesons from clover-improved lattices with Nf = 2 degen- 
erate sea-quarks - see Table HTl and Table Hill for details. 
Each of these mesons is unstable and in the Nf = 2 world 
with two degenerate quarks they have two-body decays 
to 7r7y2, TTo; and irp respectively. Here 772 is fiavour singlet, 
{uu + dd)/y2, so it is more like the ?/ than the rj meson. 
Indeed estimates .21j] of its mass from a mixture of lattice 
results and experiment suggest that it is near 0.86 GeV 
for light quarks of physical mass. Thus, for these light 
quarks with Nf — 2, the open decay channel is heavier 
for the ao meson than for the ai and 61 mesons. Hence, 
in the self-consistent world with Nf = 2 degenerate light 
quarks, we do not expect the oq meson to have any pe- 
culiar features compared to the other P-wave mesons. 

This is in contrast to quenched QCD where the flavour 
singlet pseudoscalar has the same mass as the pion, but 
an anomalous coupling. Moreover, quenched QCD al- 
lows a contribution (hairpin diagram) to the oo corre- 
lator from this two-body channel which gives significant 
unphysical effects. 

The conventional way to extract the mass of a meson 
is to use lattice simulations at successively smaller quark 
masses and to extrapolate using an expression based on 
chiral perturbation theory. For dynamical simulations, 
which are mandatory here, one has a very limited range 
of quark mass available. Resorting to partially quenched 
methods to reach lighter valence quarks is potentially 
dangerous, if the valence quarks are lighter than the sea 
quarks. Indeed a study using partially quenched methods 
on the U355 and U350 data sets has been conducted Q 
and yielded an estimate of the ao mass of 1.0(2) GeV. 
Here we explore a more reliable way to obtain the ao 
meson mass. 

Since the decay channels open to the P-wave mesons 
are quite similar, we propose to focus on the mass differ- 
ences between them since this will reduce lattice artifacts. 
The ai meson is very wide, experimentally, so that it is 
not a good point of comparison with lattice results. The 
61 meson, however is relatively narrow and should be well 
reproduced on a lattice. Indeed in ref. [ig, we were able 
to measure from the lattice the decay amplitude for the 
S-wave decay 61 — > ttuj, obtaining agreement with exper- 
iment. For lattice U355, for example, the decay threshold 
is at 0.72(4) for ao from Trri2 and at 0.883(8) for 61 from 
TTLU, in lattice units. These energy values are both above 
the mass values we report in Table iHll so each state is 
stable on our lattice and they are about equally below 
the lowest threshold. 

We show our results from two state fits to a 2 x 3 
matrix of correlators (2x2 for U350 and C390) using 
t-range 3-12 (3-10 for C390 and C410) in Table ITTTl The 
methods used are described in more detail in ref. [l8| . 
We use local and extended sources at the source (and 
two sizes of extension at the sink in some cases). The 



excited mass values are in all cases significantly higher 
(by over 50%) than the ground state values reported and 
do not correspond to any simple two-body level. Thus 
the TTri2 threshold level at aE — 0.72 for U355 does not 
feature in the fit. As we shall see later, this is consistent 
with the relatively weak transition amplitude on a lattice 
between two-body states and the ao- 

We find that the oq correlator can have big fluctua- 
tions which are apparent at large t, most noticeably for 
U350 where the zero-momentum effective mass decreases 
at large t. The origin of these fluctuations is mixing 
between the oq and the pion induced by regions of odd- 
parity in the vacuum - presumably associated with in- 
stantons. See |22l| for a discussion of oq — tt mixing in 
lattice QCD and in the instanton liquid model. With 
sufficient statistics these odd-parity fluctuations average 
to zero. Using stochastic methods (all-to-all) helps to re- 
duce these fluctuations as we reported before n^ . Using 
non-zero momentum can also act as a useful cross-check. 
This suggests that for U350 with zero momentum, we 
should use a t-range from 3-8 to reduce these fluctua- 
tion effects and retain consistency with our results from 
momentum 27rn/L where n = (1,0,0) and (1, 1,0). The 
value for U350 quoted in Table Hill is from this analysis. 

For the 6i meson, at non-zero momentum there can 
be mixing with the p meson (for some spin states). For 
the non-local (fuzzed) operators there will also be an ad- 
mixture of L = 2 (from distortion due to Lorentz boost) 
and possibly some mixing of opposite C (unless the mo- 
mentum phase factors are applied symmetrically to the 
fuzzed operator). For these reasons we rely on zero mo- 
mentum for the bi meson. 

Since the 5i meson has an unambiguous interpreta- 
tion as predominately a bound state of a quark and anti- 
quark, we show our spectrum results for it versus quark 
mass in fig ^ Here we see that our lattice results are 
quite consistent with a smooth extrapolation to the ex- 
perimental mass value for physical light quarks. To have 
a precision determination of the mass would require a 
continuum extrapolation as well as an extrpolation in 
quark mass and we do not have data sufficient to under- 
take this combined extrapolation. 

Because of the difficulties in extrapolating to light 
quarks using lattice results with a range of different lat- 
tice spacings, we focus on mass differences. Here we con- 
centrate on the difference m{bi) — m(ao) which is plotted 
against the quark mass in fig. [21 using vq determined on 
the lattices to create a dimensionless comparison. 

The point in fig. [3 labelled aq (980) assumes that the 
relevant ag meson is the lightest with mass 984.7 MeV. 
The next heaviest with mass 1474 MeV is less well es- 
tablished and would correspond to a point (-0.66) far 
below the x-axis. Our lattice results for the mass differ- 
ence show no significant dependence on the quark mass, 
and averaging our lattice results gives an estimate (using 
ro = 0.5fm) of m(6i) - m{ao) = 221(40) MeV. There 
is an additional systematic error coming from the as- 
sumption of a constant difference as the quark mass is 



Code no. 



i(7r)ro ro/a a»n(7r) 



i(p) 



C410 237 0.1410 1.29 

C390 648 0.1390 1.93 

U355 200 0.1355 1.47 

U350 151 0.1350 1.93 



3.01 0.427(1) 0.734(4) 

2.65 0.729(1) 0.969(2) 

5.04 0.292(2) 0.491(7) 

4.75 0.405(5) 0.579(8) 



TABLE II: Lattice gauge configurations U355 and U350 from 
UKQCD ^ and C390 and C410 from CP-PACS ^ are 
used, all having spatial extent L=16a. These have Nf = 2 
flavours of sea quark and we use valence quarks of the same 
mass as the sea quarks. 

Code am{bi) am{ai) am{ao) 



C410 1.17(3) 1.15(2) 1.03(4) 

C390 1.48(4) 1.39(5) 1.33(8) 

U355 0.77(2) 0.72(2) 0.64(4) 

U350 0.87(2) 0.88(2) 0.75(3) 



TABLE III: Results for P-wave mesons from the methods of 
ref. \lM for U355 and C410 and from conventional methods 
for U350 (with 4 time sources) and C390. 



decreased, which we are unable to quantify. As discussed 
above, our lattice masses in Table UTTl are at quark masses 
around the strange quark mass and at non-zero lattice 
spacing. They correspond, as expected, to masses some- 
what larger than the physical hi mass of 1230 MeV, but 
are consistent within the expected systematic errors of 
the extrapolations necessary. 

As discussed above, we do not expect the two-body 
thresholds to play a significant role in our Nf — 2 spectra. 
We do, however, measure these decay transitions to have 
a more complete analysis. 



III. HADRONIC DECAYS 

For the case of Nf — 2 degenerate quarks, the matrix 
elements for decay transitions of a non-singlet scalar me- 
son to two pseudo-scalar mesons are given in Table IIVI 
where the quark diagrams are illustrated in fig. |21 

Only one case, oq — > ^1^2, is allowed staying strictly 
within Nf = 2 with valence quarks of the same properties 
as sea quarks (here 772 is the flavour singlet pseudoscalar 
for Nf = 2, namely {uu + dd)/^/2). This case involves 
a disconnected diagram (D) and is not directly relevant 
to phenomenology. In the limit that the strange quark is 
much heavier than the u and d quarks, we expected the 
neglect of s quarks in the sea to be a good approximation. 
In that case, decays such as Kq — > Kti can be studied 
from diagram T. For the physical case with s quarks 
of some 80 MeV, oq — > ttt^s and ao -^ KK may also 
be determined adequately from Nf = 2 lattice study of 
diagram T. 

With this in mind, we first evaluate the lattice tran- 
sition amplitude corresponding to the connected trian- 
gle diagram T . The contribution of T to various decay 
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FIG. 1: Mass of the 61 mesons (in units of ro ~ 0.5fm) 
versus quark mass. The strange quark mass corresponds to 
(m(7r)ro)^ « 3.4. 



ao K K 1 

ao n rjss —1 

ao n r,s (2/3)'/" 

K* K+ Tv° 2-^/2 

K* K'' 7r+ 1 

K* K »72 2-^/2 -2-1/2 

K* K ri,s 1 -1 

TABLE IV: Coefficients of transition amphtudes from flavour 
non-singlet scalar meson S to P1P2 for the triangle quark di- 
agram (T) and the disconnected quark diagram {D). Only 
the top line is allowed if Nf — 2 strictly. The other lines 
are allowed when a valence s quark is added. We define 772 
as {uu + dd)/\^, rjss as ss and T^g as (ww + dd — 2ss)/y/6. 
We have assumed that the disconnected contributions to the 
decay to rjs cancel. 
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FIG. 2; Mass difference of 61 and ao mesons (in units of 
ro ~ 0.5fm) versus quark mass. 



amplitudes will have the numerical factors listed in Ta- 
ble ^3 The most relevant cases will be oq -^ KK and 
K* -^ Ktt. This is a partially-quenched evaluation in the 
sense that we use valence s-quarks (of the same mass as 
our u, d sea-quarks) which are not present in the sea. We 
are able to use similar methods to those used to study p 
decay [lj| and hybrid meson decay [l3| . 

The lattice results for the connected (T) contribution 
to a generic scalar meson transition to two pseudoscalar 
mesons are presented as the normalised lattice ratio 



Rit) 



(S ^ P1P2 



v/(5^5)(Pl^Pl)(P2^P2) 

where the three-point correlator is constructed from 



FIG. 3: Quark diagrams involved in the decays listed in 
Table ITVI where D is the disconnected diagram and T is the 
triangle diagram. 



propagators as illustrated for T in fig. 01 Each two and 
three-point correlator is taken at the same time separa- 
tion t. 

Since the ao mass is approximately twice the pseu- 
doscalar mass (see Table lllHl at zero momentum, we have 
an on-shell transition and we expect [13, [la| the ratio 
R{t) to be approximately linear with slope xa versus t 
where x is the lattice transition amplitude. This is in- 
deed observed, as shown in fig. fig:xtaO. 

We first checked that using different operators to create 
mesons gave essentially the same ratio R{t). We use local 
or fuzzed operators for each of the three particles involved 
and in each case the ratio is the same within errors for 
the t region of interest for the case we studied in most 
detail, namely with all momenta zero. 

The most reliable determination of the coupling con- 
stant comes from using meson operators which minimise 
excited state contributions. We use fuzzing with sepa- 
rations of 3a (C410) or 5a (U355) to achieve this. We 
extract the slope xa by taking finite differences and re- 
late this lattice transition amplitude to the continuum 
coupling via Fermi's Golden Rule. The derivation of the 
phase space factor is described in ref . [13] . Then, to com- 
pare different lattice data sets, we extract the effective 
coupling using El El 113 
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FIG. 4: The normalised ratio K(i) for the connected contribu- 
tion (T) to the transition S ^> P\P2- The contribution of T to 
particular decays can be read off from Table HVl The number 
of lattice gauge configurations analysed was 90 (U355), 165 
(C410) and 30 for each U395 case. The dotted line illustrates 
the expected behaviour with slope xa for C410. 



FIG. 5; The effective coupling g extracted from R(i) as de- 
scribed in the text for the triangle graph T for 5* — > P1P2 with 
zero momentum (also some results for non-zero momentum as 
discussed in the text). The dotted line at g = 1 is to guide 
the eye. 



Here the decay width F is, for a process with amplitude 
T, given by F/fc — g^, where k is the centre of mass 
momentum of the decay products. For particular tran- 
sitions, the quark coupling coefficients of Table Hvl also 
enter, squared, in the decay rate. 

As a first check of this approach, we evaluated the 
effective coupling from lattices that differ only in spatial 
size (labelled U395, see ref. [13 for more details) and we 
found excellent agreement when the spatial volume was 
changed by a factor of 2.4., as shown in fig. 01 

The coupling extracted, as above, from our higher 
statistics data-sets is shown in fig. O This shows a cou- 
pling g ~ 1 which has implications which we discuss later. 
The consistency between the two determinations (C410 
and U355) which have different spatial volumes and dif- 
ferent lattice spacings is satisfactory. As an overall sum- 
mary we quote a coupling g = 1.0(2). 

We also have available some results (from 40 gauge 
configurations for U355 and for 50 for C410) for transi- 
tions involving non-zero momentum, especially 5(1) — > 
Pi(0) + ^2(1) where the momentum (in units of 2tt/L) 
is given in the brackets. These results used the methods 
of ref. ]VA HsJ I respectively. The normalised lattice ratio 
R{t) is shown in fig. Eland the coupling extracted assum- 
ing the formulae above is included for C410 (where we 
used an optimum method to extract ground state con- 
tributions) in fig. 21 As discussed in ref. pj| , the decay 
in flight poses some problems of normalisation (since it 
is not quite equivalent to a centre of mass decay with 
relative momentum tt/L) , so must have a somewhat big- 
ger systematic error to compensate. Nevertheless, we see 
an approximate agreement of the lattice transition am- 
plitude xa and of the coupling g when the decay has 
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FIG. 6: The normalised ratio R{t). Here MOM=0, 1 refers to 
the transition S{k) —> Pi (0)-|-P2 (k) with momenta k = 2mi/L 
with n = Q and 1. 



momentum release of zero and of tt/L. This is to be ex- 
pected for an S-wave decay where the effective matrix 
element should be independent of momentum. 



IV. LEPTONIC DECAY CONSTANT 

The decay constants of non-singlet 0"*"+ mesons are not 
routinely calculated using lattice QCD, although they are 
of interest for a number of reasons. The value of the 
decay constant, which is basically the amplitude to find 



a quark and anti-quark at the origin, can help distinguish 
between different quark content of the meson |25ll2a |. For 
instance, if the ao was a KK molecule, then the decay 
constant would be small relative to the value of the pion 
decay constant. The decay constant of the ag meson is 
also a theoretical input to study of B meson d ecay s and of 
T decays to final states that include an aq |27ll2a . I29ll30l | . 

Diehl and Hiller discuss the prospects of determining 
the value of the decay constant of the ao mesons from 
experiment |27|. As we explain below, a direct measure- 
ment of the decay oq constant coupled with computation 
of a QCD matrix element could be used to compute the 
mass difference of the up and down quarks. 

The decay constant of the light 0+ meson can be de- 
fined by equation ^ 



(0 I V^'\ao) = ip^ga. 



(1) 



where V^ is the vector current for quark flavours a and 
b. 

The conservation of the vector current is used to relate 
the operator in equation ^ to the scalar current. 



dp-il'^lt^l'') = K'^a - mb)q''q'' 



(2) 



for light quarks with flavour a and b. This motivates a 
definition of the decay constant such as 



i{Q I q-q" I ao) = fao^l 



(3) 



To compare the size of the decay constant of the ao to 
that of the 7^*(1430) meson, Maltman ^U\ defined a new 
decay constant with a slightly different normalisation. 

The direct use of equation 12 is impossible in a lattice 
calculation with two degenerate flavours of sea quarks. 
The vector current does not couple to the scalar meson 
in this case. The decay constant in equation |21 is non- 
zero in an unquenched lattice QCD calculation with two 
flavours of sea quarks. 

There is another reason for splitting the definition of 
the 0^"*" meson decay constant into a quark mass fac- 
tor and QCD matrix element. Currently there is a dis- 
agreement between the value of the strange quark from 
unquenched lattice QCD calculations that use different 
types of fermion for the light quarks |31|. Some recent 
papers [33,|33,|3l| report summaries of the values for the 
strange quark mass published around time of the lattice 
2005 conference, using different formulations of lattice 
QCD. Lattice QCD calculations are only just starting to 
report values for the differences between the masses of 
the up and down quarks 35]. Hence, we prefer to quote 
separately our measured matrix element rather than in- 
troduce explicit factors of the quark mass. So, it is more 
natural to define the decay constant using equation^ 



(0 I qq\ao) = niagfao 
The relation between fa^ and gag is via 



"^rf 



rua 



'-f 

J a 



(4) 



(5) 



^ a fao/^S Zs 


fao MeV 


.1355 0.352(19) 0.70 


488(26) 


.1350 0.346(30) 0.71 


460(40) 


.1410 0.474(48) 0.79 


478(48) 


.1390 0.480(97) 0.84 


513(104) 



TABLE V: Our results for decay constant of the ao meson. 

The explicit factor of the quark masses {irid — rriu) in 
equation [3 is the reason that Narison 36] computes the 
value of gao to be between 1.3 and 1.6 MeV. 

The matrix element in equation 0] is extracted from 
the amplitudes in the fits to the correlators (see [23] for 
example). The raw numbers from the lattice calculation 
need renormalisation. To convert the lattice number to 
the MS scheme we use tadpole improved perturbation 
theory to one loop order ^g • The renormalisation factor 
for a scalar current, at the scale /i = 1/a, is 



Zsi^J. = 1/a) == Wo (1 - asSc) 



(6) 



where uq is the fourth root of the plaquette, and the 
constant Sc is 1.002 for the Wilson g aug e action [33 and 
0.5031 for the Iwasaki gauge action |40ll4ll]. 

To remove 0{a) terms we also need to use improve- 
ment coefficients. We define the renormalisation Zs that 
includes the improvement factor 



Zs = Zs {1 + bsrriq 



(7) 



where niq is the mass of the light quark. We used the 
one loop expression for bs- 



(1 



(8) 



where the constant bsc is 1.3722 for the Wilson ac- 
tion 113, and 1.2800 for the Iwasaki action ^,42j. We 
used the coupling computed in the MS scheme. For the 
UKQCD data set we used the coupling determined on 
the same data set [ij. For the CP-PACS data we used 
the MS coupling quoted in their paper |2J]. The cou- 
pling was evaluated at the scale /z = 1/a. Our results 
are in tabled As we only have decay constants for two 
different quark masses with the same action, we do not 
attempt a chiral extrapolation. The dependence of the 
decay constant on the pion mass seems small, however. 
In table IVll we compare our results to other determina- 
tions of the decay constants. The results in table M show 
that the decay constant fao is not suppressed relative to 
the pion decay constant. 

The decay constant of the 0++ meson is one of the 
parameters in the model that gets rid of the ghost state 
in the scalar 0"'""'" correlator in quenched QCD ^ and 
partially-quenched QCD jl^, so there are estimates for 
it. These studies of the O"*""*" used another normalisation 
convention for the scalar decay constant, so we do not 
tabulate their values here. 

Chernyak |45| uses a fit to data with a factorisation as- 
sumption to obtain 50+ = 70± 10 MeV for the AT* (1430). 



Group 


Method 


fao MeV 


Maltman [25] 

Shakin and Wang. [46j 

Narison. [361 


sum rule 

model 
sum rule 


298 

433 

320 - 390 



TABLE VI: Some results for decay constant of the ao meson. 
We used a value of nid — m„ of 4 MeV to convert the nor- 
malisation of Narison's estimate. The quark masses quoted 
by Shakin and Wang were used to convert normalisation con- 
ventions for other two results. 



Converting to our normalisation conventions, using a 
nominal value of the strange quark mass of 100 MeV, 



this corresponds to fx' 



(1430) 



1000±140MeV. There- 



suits for the decay constants in tableware larger than 
the results of UKQCD's recent calculation of the decay 
constant of the O''" charm-light meson [33 ■ 

As an aside we note that if the decay constant 170+ of 
the flQ or _R''^(1430) was measured experimentally, then it 
would allow an additional method to measure the quark 
mass differences to„ — md, iris — rud respectively, using 
lattice estimates of the QCD matrix elements. 



V. DISCUSSION 

For the non-strange flavour non-singlet scalar meson 
(flo), our lattice determinations using the self-consistent 
Nf = 2 approximation to QCD give clear support for 
a physical oq meson lying substantially lighter than the 
bi. The mass estimates we find are consistent with the 
observed ao at 950 MeV but not the heavier state at 1474 
MeV. 

To relate our approach to experiment with an addi- 
tional strange quark, we can assume that the strange 
quark pair production is relatively small and so can 
be neglected. For Kq propagation, for example, this 
amounts to treating the Ktt channel correctly but hav- 
ing an anomalous contribution from Krjss intermediate 
states. Here the ijss propagation has a missing piece (just 
as 772 does in the ?727r contribution to the oq propagation 
in quenched QCD with Nf = 2) and so will not have 
a single exponential but two contributions with masses 
corresponding to (i) the connected pseudoscalar meson 
with valence quarks of strange mass and (ii) the 772 me- 
son. Both of these contributions are not especially light, 
so we do not expect any major distortion of the Kq from 
using valence s-quarks. Similarly the ao decays to KK 
and 778''' are expected to be accessible without major dis- 
tortion from the neglect of strange quarks in the sea, as 
we discussed above. 

For the strange scalar mesons, the Kq(1430) with mass 
1412 MeV is heavier than the corresponding axial mesons 
(K* with masses 1273 and 1402 MeV). These two axial 
mesons are related to a mixture of the strange partners of 
the non-strange 61 and ai mesons since charge conjuga- 
tion is not a good quantum number for strange mesons. 



So the interpretation in this case is unclear. As well as 
this strange scalar meson at 1412 MeV, one might expect 
a lighter state, about 100-130 MeV heavier (mass split 
determined from tensor mesons) than the ao(980). The 
so-called kappa (k) at 700-900 MeV with a very broad 
width (400 MeV or more) has been claimed by many 
sources 1] and a recent analysis |47|] gives mass 750^50 
MeV. There is no consensus yet on the existence of the 
kappa, because some analyses of experimental data see 
no sign of it [ig . Our lattice studies suggest that a scalar 
Kq meson of mass around 1000-1200 MeV would be ex- 
pected in a theory with Nf = 2 sea quarks and a strange 
quark treated as a valence quark. For our case where 
the valence s-quark and u, d sea-quarks have the same 
mass, the anomalous Krjss intermediate state combines 
with the Kt: intermediate state to give only a Kri2 inter- 
mediate state, just like the case of oq propagation. Hence 
our lattice treatment does not correctly include the Ktt 
threshold in the Kq meson propagation and so may be 
less reliable than the ao propagation. 

The connected decay diagram (T) is appropriate for 
the decays ao -^ KK, qq -^ nrjs and Kq -^ Ktt where 
the appropriate factors are given in Table HVI Then the 
experimental data [l| can be used to estimate the cou- 
pling (from r/fc). For K5(1430), this giros g^ = 0.32(3). 
While for the k, one recent analysis |42| finds a width 
of 342 ± 60 MeV which corresponds^ to g^ ^ 0.7(2). For 
ao(980), the state is close to the KK threshold which 
distorts the appearance of the meson. Phenomenological 
analyses vary but one example quotes (49j a total width 
of 153 MeV and a coupHng given by g'^ = 0.82 for KK 
and around 0.7 for r77r. For ao(1450), the partial widths 
are not well known and one can only estimate that the 
KK and tjtt decays yield couplings smaller than g^ — 0.23 
and 0.34 respectively. 

Our determination of the coupling which controls de- 
cay is also relevant for identification of states. We find a 
coupling (normalised to diagram T above) given by g ~ 1. 
This favours the lighter ao and k meson over the heavier 
states. Our determination of the ao decay constant dis- 
favours a molecular structure for this state, in agreement 
with our conclusion from hadronic decays. The only con- 
cern is that for the Kq meson the experimental evidence 
gives a k meson lighter than our expectations. 



VI. CONCLUSION 

We have studied the spectrum and decay of non-singlet 
scalar mesons from first principles using lattices with 
a consistent (unitary) field theoretic interpretation for 
Nf — 2 flavours of sea-quark. Rather than extrapolate 
the scalar masses directly, we concentrate on the mass 
splitting between the ao and 61 mesons from the lattice. 
The lattice results are unambiguous and point to a scalar 
meson which is 221(40) MeV lighter than the 61. Since 
the experimental mass value of the 61 meson is 1230 MeV, 
this suggests that the aQ(980) is indeed the lightest non- 
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singlet scalar meson in a theory with Nf = 2 flavours 
of degenerate quark. Our approach does not include the 
KK channel, so this channel is to be regarded as having 
an impact on a pre-existing state, rather than as being 
the dominant component of the state. In other words, we 
do not find that a KK molecule is a good approximation 
to the ao(980). 

Our results for the decay transition amplitude are also 
consistent with the phenomenological estimates of the 
coupling of the ao(980) to KK and j/tt. Overall, we con- 
clude that the ao(980) is predominantly a conventional 
meson with normal couplings to qq. 

For the Kq scalar meson, we expect a mass 100-130 
MeV heavier than the uq (based, for example, on the 
observed mass splittings of the tensor mesons). This is 
not easily related to any experimental candidate: the k 
is too light (700-900 MeV), while the K^{U30) is too 
heavy. What may help clarification is that we find a 
decay coupling transition (to Kn) which is comparable 



to that needed phenomenologically for the k but much 
larger than that needed for the Kq{1430). This suggests 
that the k is more closely related to the state obtained 
in Nf — 2 lattice QCD. A lattice treatment with the 
strange quark included in the sea would help to clarify 
further this conclusion. 
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